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ABSTRACT 

Super star clusters (SSCs) are typically found in interacting galaxies and trace 
an extreme form of star-formation. We present a if-band study of SSC candidates in 
a sample of local luminous infrared galaxies (LIRGs) using two adaptive optics in- 
struments (VLT/NACO and Gemini/ ALTAIR/NIRI). In addition to facilitating SSC 
detections in obscured environments, this work introduces SSC studies in hosts with 
higher star- format ion rates (SFRs) than most previous studies. We find that the lumi- 
nosity functions (LFs) of the clusters are reasonably well-fitted by a single power-law 
with the values of the index a ranging between 1.5 to 2.4 with an average value of 
a w 1.9. This value appears to be less steep than the average a w 2.2 in normal spiral 
galaxies. Due to the host galaxy distances involved (median Dj_, <~ 70 Mpc) blend- 
ing effects have to be taken into account, and are investigated using Monte Carlo 
simulations of blending effects for LFs and a photometric SSC analysis of the well- 
studied Antennae system which is artificially redshifted to distances of our sample. 
While blending tends to flatten LFs our analyses show that Aa is less than <~0.1 
in our sample. The simulations also show that in the luminosity range, Mk < —13, 
considered in this work the extracted SSC luminosities are generally dominated by a 
single dominant star cluster rather than several knots of SF. We present resolution- 
and distance-dependent SSC surface density confusion limits and show how blending 
rates and aperture sizes affect the LFs. The smallest possible apertures should be used 
in crowded regions. 
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1 INTRODUCTION 

Young massive star clusters are related to triggers of star 
formation (SF) in galaxies and contain clues to the physi- 
cal conditions under which extremely strong SF happens. 
These clusters, often called "super star clusters"(SSCs) 
when masses are in the 10 5 - 10 7 M© range, are typically 
found in interacting and merging gas-rich galaxies. Their 
birth, evolution and disruption are not well understood while 
these issues are very interesting in the context of clustered 
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SF in general (e.g. Lada & Lada 2003) and in particular since 
they might be the progenitors of globular clusters (e.g. Ash- 
man & Zepf 1992; Holtzman et al. 1992; Ho & Filippenko 
1996; Elmegrecn & Efrcmov 1997). For a recent review of 
young and massive star clusters see Portegies Zwart et al. 
(2010). 

SSCs can be used to trace the history of bursts of SF 
in galaxies (e.g Bastian 2008; Escala & Larson 2008; Escala 
2011; Adamo et al. 2011b) since SSCs are bright and rel- 
atively simple to model as single stellar populations. It is 
important to understand how wide-spread starbursts can be 
in interactions and mergers, particularly in the context of 
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high-rcdshift galaxy formation, since these conditions could 
be similar to what is happening at z>l. Local LIRGs (lumi- 
nous IR galaxies, logL/ij/L© > 11) exhibit extreme star for- 
mation and large populations of SSCs and they are believed 
to be good analogs for higher-z star formation in general 
(e.g. Alonso-Herrero et al. 2009; Elbaz et al. 2011; Tekola 
et al. 2012). 

Although the young and compact SSCs (age ~ 10- 
100 Myr, r eff ~3-5pc) are often confined to starbursts and 
interacting systems (e.g. Whitmore et al. 1999) they can 
also form in more quiescent environments such as in the 
circum-nuclear star-forming rings (e.g. Benedict et al. 1993), 
in nearby dwarf galaxies (Ostlin et al. 1998) and even in nor- 
mal spiral galaxies (Larsen 2002). The difference is just that 
the population of SSCs is much larger in galaxies exhibiting 
an extreme starburst environment with a violent SF. Over 
the past decades, research has focused more on SSCs within 
nearby (Dl < 25 Mpc) non-LIRG starburst systems, simply 
because it is much easier (e.g. Whitmore et al. 1993, 1999; 
Bik et al. 2003; Haas et al. 2008). SSCs were, in fact, one 
of the first discoveries of the Hubble Space Telescope (HST) 
while imaging the central regions of NGC 1275 (Holtzman 
et al. 1992). In particular the closest major merger system, 
the Antennae (NGC 4038/4039, log L IR = 11.0) has been 
amongst the most studied SSC hosts, providing a sample 
of thousands of these clusters. There have been some in- 
dividual cases studied further away, and at higher SF lev- 
els, such as Haroll with logL/jj = 11.22 at 81.2 Mpc and 
the Bird galaxy (IRAS 19115-2124) with logL /fl = 11.87 at 
206 Mpc, by Adamo et al. (2010) and Vaisanen et al. (2008), 
respectively. The only published work thus far characterising 
SSCs from a significant sample of galaxies is that of Miralles- 
Caballero et al. (2011), who analysed optically-selected star 
forming knots, SSCs, or complexes of clusters in 32 LIRGs 
and ultraluminous IR galaxies (ULIRGs, log Lir/Lq > 12). 
Their results indicate for example that SSCs in post-mergers 
are larger, more luminous and redder, than in less advanced 
systems. An even larger sample of SSCs in 87 LIRGs has 
been recently studied by Vavilkin (2011). 

It is extremely important to study SSC characteristics 
over a wide range of galaxy properties and SF environments, 
given that there is a very lively debate on the universality 
of SSC mass functions and disruption mechanisms (see e.g. 
Lamers 2009). For example, Fall et al. (2009), Whitmore 
et al. (2007, 2010) and Chandar et al. (2010, 2011) argue 
that any turnover in the mass-function results from a con- 
stant disruption triggered by internal mechanisms (e.g. stel- 
lar evolution, 2-body relaxation) over the cluster lifetime. In 
contrast, Lamers et al. (2005), Gieles et al. (2006c), Larsen 
(2008), Bastian et al. (2011) and Kruijssen et al. (2012) pro- 
pose that the mass of the cluster and/or the environment 
of its host galaxy play an important role in the disruption 
model (low-mass clusters are less likely to survive the dis- 
ruption, while environments with low giant molecular cloud 
densities will endure less disruption) that will later define 
the cluster mass function (CMF) shape. 

The luminosity function (LF) of SSCs, which is sug- 
gested to be a reflection of the cluster initial mass function, is 
a common tool used to help understand their formation and 
evolution (e.g. Whitmore et al. 1999; Gieles et al. 2006a). 
While theoretical and observational studies have agreed that 
the LF of globular clusters is well-fitted with a Gaussian 



function, there are varying results regarding the LF of SSCs: 
though power-laws are usually found, and expected theo- 
retically, log-normal distributions and bent power-laws have 
also been advocated (e.g. Elmegreen & Efremov 1997, Fritze- 
v. Alvensleben 1999, Larsen 2002, de Grijs et al. 2003, Lada 
& Lada 2003, Gieles et al. 2006a, Anders et al. 2007, Whit- 
more et al. 2010, Miralles-Caballero et al. 2011). To compli- 
cate matters, there are several observational issues to over- 
come as well to pin down the shapes, relating, for exam- 
ple, to stellar contamination, completeness, and blending 
(e.g. Robin et al. 2003; Anders et al. 2007; Fall et al. 2009; 
Adamo et al. 2010). The typical single power-law approxi- 
mation (e.g. Elmegreen & Efremov 1997) takes the form: 

N(L)dL ~ L~ a dL, (1) 

where a~2 with scatter in the range 1.5 to 2.5. Gieles 
et al. (2006b), Bastian (2008), Haas et al. (2008) and Larsen 
(2009) among others, on the other hand, suggest a double 
power-law or a Schechter function to better fit the cluster lu- 
minosity function (CLF) due to a turnover located at higher 
luminosities, which could be a sign of a truncated CMF. 
The interesting question is whether there are any system- 
atic trends to those LF index values. Normal spiral galaxies 
tend to have LFs biased to the steeper range of slopes at 
a ~ 2 - 2.4 (e.g. Larsen 2002; Gieles et al. 2006a; Mora 
et al. 2009). Recent results from larger samples of LIRGs 
(Miralles-Caballero et al. 2011; Vavilkin 2011) show evidence 
that a w 1.8 — 2.0, though at distances > 200 Mpc where res- 
olution and blending effects may play a role, they find that 
the LF flattens toa~ 1.7. Moreover, Ostlin et al. (2003) and 
Adamo et al. (2010, 2011a,b) derived flatter slopes rang- 
ing between 1.5 to 1.8 for the /-band CLF in their very 
strongly star forming blue compact dwarf galaxies (54 to 
82 Mpc), while Surace et al. (1998) and Inami et al. (2010) 
find similar slopes between 1.1 and 1.8 in the B-band for 
three (U)LIRGs. 

In our ongoing study we are expanding the SFR range 
of SSC host galaxies to include LIRGs from their lowest 
limit to ULIRGs, and including galaxies in various stages of 
interactions, from isolated and paired galaxies, to interact- 
ing, merging and merger remnant stages with a sample of 
(ultimately) dozens of targets (Vaisanen et al. 2012). Fur- 
thermore, most of the studies thus far have been done in 
the optical, making extinction effects potentially difficult in 
the notoriously complex dusty environments of gas-rich in- 
teractions. We use iCs-band 1 adaptive optics (AO) observa- 
tions which match the spatial resolution of the HST optical 
studies. Crucially, the use of near infrared (NIR) has great 
potential in opening a new angle into the SSC populations 
in that they probe deeper into the dusty birth regions of 
the SSCs and the very obscured regions in the inner parts 
of the galaxies. There is also an interesting time- window at 
ages of ~ 10 Myr when SSCs are expected to be very NIR- 
luminous due to their high-mass stars entering the red su- 
pcrgiant (RSG) phase. There are very few NIR studies of 
SSCs compared with the optical ones, including Lai et al. 
(1999), Mengel et al. (2005), Alonso-Herrero et al. (2006), 



1 Hereafter, we will refer to both the Johnson K-band and 
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Pollack et al. (2007), Vaisanen et al. (2008), and Adamo 
et al. (2010). 

The main specific objectives of this work are to derive 
the Jf-band LFs of the massive star clusters and to evaluate 
the effect of blending on the LFs. In forthcoming papers we 
will study the mass functions and age distributions of the 
targets by combining optical data with the NIR photometry. 
The paper is organised as follows: we describe the data and 
its reduction in § 2. The cluster analysis is presented in § 3. 
The NIR CLFs of the sample are presented in §4 and are 
interpreted in § 5; finally, we summarise our findings then 
suggest our future work in § 6. 

Throughout the paper we assume the standard cosmol- 
ogy: H = 73 km s" 1 Mpc -1 , Q M = 0.27, and ft A = 0.73. 



2 OBSERVATIONS AND DATA REDUCTION 

We analyse a representative sample of ten local LIRGs from 
an ongoing NIR AO survey which is mainly intended to 
search for dust-obscured core-collapse supernovae (see e.g. 
Mattila et al. 2007; Kankare et al. 2008, 2012; Vaisanen 
et al. 2009). The targets analysed here were drawn from 
the IRAS Revised Bright Galaxy Sample (Sanders et al. 
2003). They lie at redshifts of 0.OKK0.O5, their lumi- 
nosities range from \og{LiR/Lo) = 11.2 to 11.9, and they 
were selected to be starburst-dominated based on their IRAS 
colours (/25//60 < 0.2). The observations were taken in K- 
band using two different ground-based instruments with AO 
imaging: the NAOS-CONICA on the ESO Very Large Tele- 
scope (VLT/NACO) and the ALTAIR/NIRI on the Gemini- 
North telescope (ALTAIR/NIRI). For AO correction, natu- 
ral guide stars were used for the NACO data, whereas laser 
guide stars, with a tip/tilt reference star, were used with the 
ALTAIR/NIRI data - this requirement of suitably bright 
reference stars is a further sample selection constraint, but 
does not bias the LIRG characteristics in any way. 



2.1 NACO data 

Our present sample includes two NACO galaxies, 
IRAS 18293-3413 and IRAS 19115-2124, the latter also be- 
ing the subject of our pilot SSC study in relatively distant 
galaxies compared to most SSC works (Vaisanen et al. 2008). 
The data taken with the VLT UT4 and NACO S27 cam- 
era have a plate scale of 0.027 "pix -1 and a field of view 
(FOV) of 27". The AO correction worked well resulting in a 
FWHM ~ 0.1" for the point sources. Frames were taken with 
exposure times of 30 sec in dithering mode with an integra- 
tion time per pointing of 90 sec, and Table 1 lists the total 
integration times. More details of the NACO observations, 
as well as data reduction, are given in Mattila et al. (2007) 
and Vaisanen et al. (2008) for the two galaxies, respectively. 



2.2 ALTAIR/NIRI data 

The other 8 targets come from a recent multi-epoch sur- 
vey using Gemini-North during 2008-2012. The pixel scale 
is 0.022"pix~ 1 , yielding a FOV of 22". These data also have 
a final resolution of ~0.1". Each individual frame has an 



exposure time of ~ 30 sec, and in this case separate sky- 
frames were taken for sky subtraction in addition to 011- 
target dithering. Refer to Table 1 for details of observations. 

Our Gemini data were reduced using IRAF 2 -based tasks 
including flat-fielding and sky subtraction. Individual frames 
with significantly lower quality PSFs were excluded, as well 
as some with abnormal electronic noise. A weaker, horizon- 
tal stripe pattern was still evident and was removed by a 
custom-made de-striping algorithm. The final images were 
produced by co-adding by average-combining the individual 
frames from different observing runs after shifting them to 
a common reference. Table 1 again lists the effective integra- 
tion times per target. The astrometry calibration was per- 
formed with the IRAF task CCMAP: we downloaded archival 
H ST /ACS data of the fields and accurately re-calibrated 
these using the Guide Star Catalog II, and then added the 
World Coordinate System (WCS) into the FITS headers of 
our NIR images using the larger FOV HST frames as refer- 
ence images. 

2.3 NOTCam data 

We obtained if-band images of all the Gemini targets using 
the Nordic Optical Telescope (NOT) NIR Camera and spec- 
trograph (NOTCam) . In the wide field imaging mode, NOT- 
Cam has a pixel scale of 0.234"pix _1 and a FOV of 4'. The 
images were flat-field corrected, sky subtracted and com- 
bined using IRAF-based tasks. The NOTCam images were 
used as intermediate images to determine photometric zero- 
points for the smaller FOV NIRI images. 



3 THE SSC CANDIDATES 

3.1 Object detection and photometry 

For object detection we ran SExtractor v2.5.0 (Bertin 
& Arnouts 1996) on unsharp-masked versions of the im- 
ages (see the top-right panel of Figure 1 in the case of 
IRAS 18293-3413); unsharp masking was done in order to 
make point-source detection more uniform in varying back- 
ground conditions, though we stress that photometry (see 
below) was done always on original images. Critical param- 
eters for detection were tuned to minimise spurious sources 
and to include faint and also extended sources in the out- 
put catalogues. A threshold of ~ 1.5a above the background 
RMS noise combined with an upper value of the minimum 
number of ~ 10 adjacent pixels above threshold were even- 
tually chosen. 

We performed aperture photometry on the combined 
images using the task IRAF/PHOT with 3 and 5 pixel aper- 
ture radii (approximately equal to a FWHM ~ 0.1" of a 
point source) and sky annuli from 5 to 7 and 7 to 10 pixels 
(~ 0.06" wide) for the NACO and the ALTAIR/NIRI data, 
respectively. A small sky annulus is necessary for a good 
sky sampling in the strongly varying background in between 
larger-scale features of the galaxy, while small apertures are 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Galaxy name 


Exp time 
(sec) 


RA 
(J2000) 


DEC 
(J2000) 


I 

(degrees) 


b 

(degrees) 


log L IR 

(£©) 


SFR 

(MQyr- 1 ) 


rn - M 


D L 
(Mpc) 


IC 694 


1260 


11 28 33.5 


+58 33 45 


141.9 


55.4 


11.60° 


67.7 


33.28 


45.3 


NGC 3690 


2192 


11 28 29.8 


+58 33 43 


141.9 


55.4 


11.52" 


56.3 


33.28 


45.3 


CGCG 049-057 


1680 


15 13 13.1 


+07 13 32 


8.9 


50.9 


11.27 


31.6 


33.76 


56.4 


TRAS F1 7578-0400 


1470 


1 8 00 31 9 

1U UU O -L ■ U 


— 04 00 53 

VJ L ± \J\J at) 


23.4 


9.4 


11.35 


38.1 


33.79 


57.3 


IRAS F17138-1017 


990 


17 16 35.8 


-10 20 39 


12.2 


15.6 


11.42 


44.7 


34.29 


72.2 


IRAS 18293-3413 


1230 


18 32 41.1 


-34 11 27 


0.1 


-11.3 


11.81 


109.8 


34.37 


74.6 


MGC +08-11-002 


1140 


05 40 43.7 


+49 41 41 


161.6 


9.9 


11.41 


43.7 


34.51 


79.9 


IR AS F16516-0948 


900 


16 54 24.0 


-09 53 21 


9.5 


20.5 


11.24 


29.5 


34.88 


94.8 


IC 883 


1440 


13 20 35.3 


+34 08 22 


82.9 


80.6 


11.67 


79.5 


35.02 


101 


IRAS 19115-2124 


1410 


19 14 30.9 


-21 19 07 


16.1 


-14.4 


11.87 


126 


36.56 


206 



Table 1. Total exposure times, equatorial (RA, DEC) and galactic (I, b) coordinates of the observed sample are shown in this table 
ordered by distance. The luminosity distance, Dt,, and the distance modulus, m — M, are retrieved from NASA/IPAC Extragalactic 
Database (NED) while the logarithmic value of the IR luminosity, Ljr, is as estimated by Sanders et al. (2003) with a slightly different 
setting of the cosmological parameters. a These two galaxies are components of an interacting system Arp299; the total Ltr is 

divided between IC 694 (known as nucleus A) and NGC 3690 (nuclei B+C) following Mattila et al. (2012) the circum-nuclear SF portion 
being evenly split. An empirical relation using Ltr by Kennicutt (1998) is used to derive the SFR. 




Figure 1. The NACO field of IRAS 18293-3413. Top left: The field of the whole interacting system with contours demarcating the 
four selected background regions for completeness analysis. Top-right: A slightly smaller field around the primary galaxy is shown after 
unsharp masking. The SSC detections are made from this image. Lower left: All SSC candidate detections are overlaid as white on the 
original image (where the photometry is performed) while those that meet all the SSC candidate selection criteria presented in Section 3.2 
are shown as red. Lower-right: Same as previous, but only for a zoomed-in region around the nucleus. Small tick marks are in 1" units, 
except in the last panel where they are in 0.1" units. 
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needed to minimise blending effects, especially in the case of 
the crowded SSC populations detected in the NACO data of 
IRAS 18293-3413 (hence the smaller apertures with NACO 
data) . 

Since the apertures are small, aperture corrections are 
essential. In addition, the PSF shape is expected to vary 
across the frames as a function of distance from the star used 
as the AO-reference. Aperture corrections were determined 
based on the curve-of-growth (out to 1.0") of sufficiently 
bright and isolated stars at various locations throughout the 
images. In the case of NACO data the aperture corrections 
were found to be dependent on the distance of the AO ref- 
erence star, while in the Gemini data, where the laser guide 
star is located in the centre of the field within the target 
galaxy, no such systematics were found. These aperture cor- 
rections a c , as a function of radial position from the AO 
star if required, were then applied to measurements of can- 
didate clusters using an aperture of 0.1" radius. The NACO 
data a c values applied to individual candidate SSCs range 
from -1.23 to -2.07 mag and for the ALTAIR/NIRI data 
we adopted an averaged constant a c = —1.23 mag. We esti- 
mate the uncertainty of the aperture correction in a single 
frame to be typically ~ 0.3 mag. Note that these aperture 
correction methods implicitly assume that we are generally 
detecting point sources, which is supported by our simula- 
tions described in Sections 3.2, 3.3 and5.1. 

Finally, photometric zero-points were referred to the 
2MASS /fg-band point source catalogue using 2MASS com- 
mon stars in the FOV, making our magnitude system the 
same as that of the 2MASS Vega-based Ks filter. In cases 
where the 2MASS point-sources were outside the small FOV 
of the Gemini data, NOTCam if -band images were used 
as intermediate images to estimate the zero-points. For 
MCG +08-11-002 and CGCG 049-057, comparing the total 
integrated flux of the galaxy with the corresponding 2MASS 
magnitude was the only option to estimate their zero-points. 
We estimate the absolute calibration to be accurate to ~ 0.4 
mag including 2MASS catalogue photometric uncertainties 
and systematic errors in our photometry and aperture cor- 
rections. 



3.2 Cluster candidate selection 

Deciding whether a detected source is a SSC candidate or 
not is challenging when working with a single filter and deal- 
ing with a sample of distant targets where both stars and 
star clusters have similar PSF sizes. In our case, the follow- 
ing steps were carried out to generate the final if -band SSC 
catalogue for each target. 

(i) First of all, only objects falling on detectable optical 
or NIR emission from the galaxy within the frames were 
considered. The edges of frames where the noise is higher 
were also excluded in case the host galaxy extended there. 
After this process, we assumed that the catalogue is mainly 
composed of real sources since we had already tuned the 
detection parameters to avoid false detections in the inner 
regions of the frames. 

(ii) Foreground contamination by Milky Way stars is 
a real possibility, especially for targets with low galactic 
latitude (6<|20°|, see Table 1). Therefore, we estimated 
the potential effect of contamination in our data by using 



Name 


N(SSC) 
(<r m «S 0.25) 


Comp. limit (JsT-band) 
App mag Abs mag 


IC 694 


81 


19.7 


-13.6 


NGC 3690 


81 


19.4 


-13.9 


CGCG 049-057 


10 


19.4 


-14.4 


TEAS F1 7578-0400 


45 


20.4 


-13.4 


IRAS F17138-1017 


60 


20.0 


-14.3 


IRAS 18293-3413 


204 


19.9 


-14.5 


MCG +08-11-002 


12 


19.2 


-15.3 


IRAS F16516-0948 


41 


19.9 


-14.9 


IC 883 


29 


20.2 


-14.8 


IRAS 19115-2124 


56 


20.1 


-16.5 



Table 2. After imposing our selection criteria, the final number 
of SSC candidates for each target is given. We tabulate also the 
apparent and absolute magnitudes of the 50 % completeness level 
corresponding to the background region where more than m 50 % 
of the data points are below the contour level limiting that region 
(see text). The brighter absolute magnitude completeness limit 
for IRAS 19115-2124 is due to its much larger distance compared 
to the other targets. 



the Besangon model (Robin et al. 2003). IRAS 19115-2124, 
IRAS 18293-3413, and IRAS 17578-0400 yield the largest 
numbers of expected if -band stars according to the model, 
but even in these fields the maximum numbers of contam- 
inating stars at K ~ 16 mag corresponding to the brightest 
magnitude bins of SSCs, are at most at ~ 10 % of the total 
SSC candidate numbers within a typical 20" x 20"area of 
the galaxies. At fainter magnitudes star counts are at negli- 
gible levels, as are counts at all magnitude bins for the rest 
of the fields, according to the models. However, as a consis- 
tency check, we also counted obvious point sources in the 
sky areas outside of the galaxies in our fields, derived their 
surface densities in magnitude bins, and compared to ob- 
ject numbers within the host galaxy area. As an example we 
show the low galactic latitude IRAS 18293-3413 case in the 
left panel of Fig. 2. These numbers are consistent with the 
Besangon model predictions at bright magnitudes, but tend 
to be somewhat larger at K > 16: it could, for example, be 
that there are some real SSCs well outside of the projected 
galaxy area of our targets, or that the models do not probe 
Milky Way stars at faint enough levels. Though there might 
be some contamination at the brightest K ~ 15 — 16 bins, 
these would be "obvious stars" that are removed in the fi- 
nal selection step (v) below and at fainter levels the star 
numbers still are not large enough to make any significant 
difference to the ultimately derived luminosity functions of 
SSCs. We hence conclude that foreground contamination is 
not significant for our results, and no further star vs. SSC 
separation was attempted for our dataset. 

(iii) Photometric uncertainties were taken into account 
fairly conservatively to have a robust list of SSC candi- 
dates: we excluded all detections having errors <r m > 0.25 
mag. We note that excluding less secure detections can in 
principle bias eventual LF values since fainter objects will 
not be included. However, we checked that if we had selected 
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Figure 2. Left: Raw surface densities of detected objects within the galaxy area of frame IRAS 18293-3413 are shown in blue, while 
those of objects found outside the galaxy area in the same frame are shown in red. As the latter are likely foreground stars, it is seen 
that stellar foreground contamination is insignificant even in this low galactic latitude case. Apparent A"-band magnitude scale is shown 
at the top. Right: An example of SSC candidate selection for the same field of IRAS 18293-3413. A FHWM vs concentration index C 
plot of the detected objects from SExtractor with the cutoffs for candidate selection indicated. The limits are checked with simulations 
to be consistent with the parameter space of extracted intrinsic point sources (~ 4.5 pix in this case) in the fields. The filled black points 
represent the final SSC candidates, whereas all detected objects are shown in green. 



a m < 1.0 mag sources instead, we would have had only be- 
tween 6 to 14% more SSC candidates, and the final LF 
shapes and slope changes would have been negligible within 
the completeness levels considered in the analysis. 

(iv) We excluded objects that appeared to be extended 
in order to remove background galaxies and resolved HII 
regions, while also taking into account the fact that fainter 
point sources will appear to have a broader FWHM when 
extracted on top of the varying diffuse background of the 
galaxy. For a candidate to be included in our photometric 
catalogue, we imposed cuts on the values of a "concentration 
index" C, defined as the difference in magnitudes between 
the 3 and 5 pixel radius apertures (m,3 px — m$ px ), as well as 
cuts on the FWHM. A small value of the index (C ~ 0.5) 
typically results from a high S/N point source, such as a 
foreground star, whereas larger values (C ~ 1.0 and above) 
indicate potentially extended objects. Since the value of C is 
less accurate for fainter sources, we also used the FWHM of 
the detection, as measured by the task RADPROF in pixels to 
make the selections more robust. Plots such as that shown 
in the right panel of Figure 2, were generated for each target 
field to help us distinguish between truly fuzzy objects and 
unresolved SSC candidates at different S/N levels. 

The cutoff values were decided with the help of simula- 
tions and also by-eye checks interactively for consistency. As 
part of Monte Carlo simulations to define completeness cor- 
rections (described in Section 3.3) simulated intrinsic point 
source PSFs were extracted from the real data frames, and 
the cutoffs in FWHM and C were adjusted to encompass the 
output parameter space of the detections of the input point 
sources in the simulations made on the real data frames. 
The differing AO corrections on different datasets introduces 
some variation, but in most of the cases SSC candidate se- 
lection included objects having values of 0.5 < C < 1.0 and 
0.6 < FWHM < 10.2 pixels. In the case of IRAS F16516- 
0948 all objects are generally elongated, possibly due to 
a degraded AO-correction; however, since all analysis was 



done using PSFs constructed from bright stars on the frame 
itself, this does not bias the selection in any way. 

Note that retaining only unresolved point sources at this 
step does not exclude the possibility of blending discrete 
SSCs smaller than our resolution PSF into one detected SSC 
candidate - this is discussed in more detail in Section 5.1. 

(v) Finally, we made a visual inspection of the images 
with the selected candidates in order to remove the nucleus 
of the galaxy from the final catalogue and to remove several 
obvious very bright stars which happened to fall on top of 
the galaxies. 

After applying the selection criteria detailed above we 
obtained the final catalogue for each target field. The total 
number of SSC candidates per LIRG are listed in the sec- 
ond column of Table 2. Figure 1 shows, in the lower panels, 
the spatial distribution of the SSC candidates in the case 
of IRAS 18293-3413: all the sources detected originally by 
SExtractor are marked on the image, but only those which 
met all the selection criteria above are in red (most are ex- 
cluded because of their magnitude uncertainty in this case) . 



3.3 Completeness correction 

The raw SSC catalogues obtained in the previous section 
are affected by incompleteness due to photometric detec- 
tion limits of the observations; the varying and complicated 
background resulting from the different diffuse components 
of the host galaxies; and potential crowding of sources. To 
correct the data for incompleteness bias, we ran Monte Carlo 
(MC) completeness simulations with each science image. A 
PSF model for each target field was first created. The LIRGs 
MCG+08-11-002 and CGCG 049-057 did not have isolated 
stars in their fields, and for them we used a representative 
PSF model from other fields having a similar distance from 
their tip-tilt reference star, and hence similar AO-correction 
based from the reference star. Using the model, we created 
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Figure 3. The results of Monte Carlo completeness simulations 
for IRAS 18293-3413 within regions of different background levels. 
The green dashed line corresponds to the innermost region with 
the highest background. Upper panel: The fraction of simulated 
point sources recovered as SSC candidates in the simulation as 
a function of apparent and absolute K magnitude in the lower 
and upper axes, respectively. The 50 % and 80 % completeness 
limits corresponding to the "middle2" region are a reasonable 
approximation of overall completeness limits of point sources in 
IRAS 18293-3413; they are shown as the horizontal dashed lines. 
Lower panel: The y-axis plots the input minus output magnitude 
as a function of the input magnitude in the same MC simulation. 
The error bars reflect the scatter of this difference in the simu- 
lation and is plotted for one curve only for clarity. The vertical 
dashed lines show the two completeness levels determined above 
from the "middle2" region. 



artificial stars with IRAF/DAOPHOT which were used in the 
simulation. 

The main idea of the simulation is to record how many 
of the input objects, added randomly to the science images, 
are detected by SExtractor at a given magnitude range with 
exactly the same manner and configuration parameters as 
the ones chosen for optimal detection in §3.1, as well as 
with identical selection criteria to be included as an SSC 
candidate. As mentioned above, these simulations were also 
used to define the latter selection criteria since the input 
source is by definition a point-source. We also checked the 



simulations for varying PSF sizes and verified that PSFs in 
a realistic range of point source sizes in the target fields do 
not result in corrections which would change any of the main 
results concerning the shapes of the eventual LFs. 

The simulation ran from 13 to 22 if-mag in steps of 0.5 
mag. We generated 1000 random positions of artificial star 
centroids. To avoid systematic errors, we used a new subset 
of random positions for each magnitude step. In addition we 
ran the simulation separately at different background levels 
of the field. These levels were determined by defining three 
(four in cases of IRAS 19115-2124 and IRAS 18293-3413, the 
latter having the largest dynamical range between the back- 
ground and point source brightness) approximately equal 
ranges in the pixel values of a smoothed background map in 
a logarithmic scale, ranging from an essentially empty sky to 
the LIRG's core, but excluding the nucleus itself. Figure 3, 
upper panel, shows as an example the resulting complete- 
ness curves for the case of IRAS 18293-3413. The simulations 
were also used to test the accuracy of the photometry for 
systematic effects. The lower panel of Fig. 3 shows the differ- 
ence of input and output magnitudes of the detected sources. 
In the fainter bins, magnitudes cannot be measured reliably 
anymore due to very variable background, crowding, and 
inaccurate aperture corrections. However, the Amag values 
become larger than our typical overall photometric uncer- 
tainties only at or below the completeness limits relevant to 
each galaxy and background region used for subsequent LF 
analysis, and any corrections for these possible systematica 
are not attempted. 

Each individual SSC photometric data point was then 
corrected for the incompleteness bias with respect to its ob- 
served magnitude as well as its location in the complex back- 
ground field, as long as its magnitude value was above the 
50 % completeness limit of the region in question. Represen- 
tative completeness values are shown in Table 2. After this 
correction, we are ready to construct the luminosity func- 
tions of the SSCs in each target field. 



4 CONSTRUCTING THE A'-BAND SSC 
LUMINOSITY FUNCTION 

Once the SSC candidates have been selected for each tar- 
get, we can construct a binned LF and then fit a function 
to its shape. LFs are constructed as a function of absolute 
Mk magnitude; the distance modulus m — M of the ob- 
served sample listed in Table 1 was used in the conversion. 
Completeness corrections were applied to the counts of SSC 
candidates as a function of observed magnitude and back- 
ground region of the detections, as discussed above. Figure 4 
shows all the observed LFs of our targets; both raw and 
completeness-corrected counts are shown. 

After this, we fit a power-law of the form N(L)dL ~ 
L~ a dL in log-log space to the LF shape. The fitted data 
points are weighted using their respective Poisson-noise un- 
certainties y/N. In a mag-logN plot, a linear fit is expressed 
as follows: 



logiV(Mjf) = P M K + con 



(2) 



where the relation between the power-law index a (see Eq. 1) 
and the linear slope f) is (Elmegreen & Efremov 1997): 



a = 2.5/3 + 1 



(3) 
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Figure 4. ftf-band LFs of the sample using a constant bin size. The black solid line is the incompleteness-corrected LF while the dashed 
grey line is the original one. The single power-law fit of the data points is represented by the cyan solid line while, in the case of 
IRAS 18293-3413, the grey solid one and the dashed blue line result from broken power-law fits, and the curved green line is a Schcchtcr 
function fit. The vertical lines mark the 50% (grey) and the 80% (cyan) completeness levels. The y-axes scales are not the same, because 
the number of SSC candidates differs from one target to another, making the optimal constant bin size different in each case. 
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Name 


a con 


a var 


X 2 , 

A - r e a 


IC 694 


2.29±0.09 


2.36±0.10 


0.85 


NOC 3690 




i so+n in 


0.43 


IRAS F17578-0400 


1.74±0.11 


1.78±0.07 


0.82 


IRAS F17138-1017 


1.53±0.12 


1.64±0.08 


0.58 


IRAS 18293-3413 


1.88±0.15 


1.89±0.09 


1.49 


IRAS F16516-0948 


2.17±0.15 


2.41±0.07 


2.17 


IC 883 


1.45±0.13 


1.58±0.07 


0.47 


IRAS 19115-2124 


2.08±0.10 


1.97±0.08 


0.61 


Average: 


1.87±0.30 


1.93±0.31 





Table 3. Power-law indices from weighted linear fitting of the 
LFs in Fig. 4 using the relation in Eq. 3. a\ on and a" ar are, re- 
spectively, the indices derived from binning with a constant and a 
variable bin width. Xred show the reduced Chi Square values for 
the single power-law fits using the constant binning. The uncer- 
tainties in the slopes a are derived from the rules of propagation 
of errors in Eq. 3, after calculating the uncertainty in /3 which is 
the weighted linear slope shown in Eq. 2. Note that in the case of 
a small dataset, e.g. IRAS F16516-0948 and IC883, the value of 



may not be a good representation of the goodness of the fit. 
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Figure 5. A"-band SSC LF of the sample, except the dataset from 
IRAS 19115-2124. The single power-law fit of the data points is 
represented by the cyan solid line while the vertical lines mark 
the average values of the 50 % (grey) and the 80 % (cyan) com- 
pleteness levels. 



Due to the very small number of SSC candidates in 
MCG+08-11-002 and CGCG 049-057, N=12 and 10 respec- 
tively, we did not fit their LFs. 



4.1 Single power- law fits 

As a first step we fitted the SSC candidate LF distribu- 
tions with a single power-law shape using a constant bin 
size. Since the targets do not have the same number of SSC 
candidates, each galaxy has its own constant bin size to bal- 
ance having statistically enough sources per bin while try- 
ing to maximise the number of bins overall. The LF bins 
were fitted from the brightest bin down to the last bin 
above the 80% completeness limit, the data points plot- 
ted as squares with error bars in Fig. 4 indicating the bins 
which were fit. The resulting power-law indices estimated 
using Eq. 3 are shown in Table 3 as af n , and the fitted 
LFs are plotted in Fig. 4 as the cyan line. The values of 
af™ range from 1.5 to 2.4. The average over the sample is 
a\ on = 1.87±0.30, or af" = 1.84±0.28 if the most distant 
LIRG, IRAS 19115-2124, is excluded. The quoted errors are 
the formal uncertainty of the fits. Fitting a single power- 
law function to the LF of the combined dataset down to a 
Mk — —14.5 mag (which includes the SSCs from all 8 LIRGs 
in Fig. 4 except for the most distant target IRAS 19115-2124 
which has a significantly different magnitude), gives a slope 
a con = 1 g2 ± o 25 (pig. 5). This is quite similar to the av- 
erage slope. 

To check whether a single power-law is a good approx- 
imation of the LF, we estimated the reduced chi-square 
statistic Xred (the ratio of chi-square x? an d the degrees 
of freedom of the dataset) in each case. In most cases Xred 
values (see Table 3) indicate that a single power-law appears 
to be a reasonable fit to the data. Note that the value of Xred 
may not reflect the goodness of the fit when dealing with a 
small number of data points. 



4.2 Effect of binning on the LFs 

The shape of the LF may be affected by sample binning: 
for example, Mafz Apellaniz & Ubeda (2005) argued that 
the difference in the value of a can be as large as 0.3 for 
small datasets. To check for consistency of our results, we 
constructed the LFs using two different methods: first using 
a constant bin size as above, and secondly using a variable 
bin size and assigning an equal number of objects to each bin, 
as proposed by Mafz Apellaniz & Ubeda (2005). Table 3 lists 
the values from both methods: a\ on is the slope of the LF 
for a constant bin size, while a\ ar is the one that results 
from using a variable bin size. Though the values derived 
from the two binning methods are slightly different, they 
are consistent within the error estimates of each fit. 

In addition, we ran systematic tests on the effect of 
using a particular size of the constant-size bin. The scatter 
of the fitted index values a with a wide range of bin sizes 
is of the same order, ~ 0.15, as the uncertainties of the 
slope fits presented in Table 3, and we conclude that bin 
size does not influence our final results. The characteristics 
of the resulting LF slopes will be discussed in § 5. 



4.3 IRAS 18293-3413: Schechter and broken 
power-law fits of the CLF 

Figure 4 suggests that in the case of IRAS 18293-3413 a two- 
component power-law would be a better fit to the CLF. 
Hence we also fitted this LF of the target with a broken 
power-law with both slopes as free parameters, but the 
break, or bend point, Mbend fixed at a location chosen by 
eye. The form of the fit is expressed as follows in a mag-logN 
plot: 



logN(M K ) 



/3 3 M K 
f3 2 M K 



for 
for 



M K > M bend 
M K < M bend 



(4) 
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where fiz and fiz are related, respectively, to the bright and 
faint slopes a 2 = 2.17±0.15 and a 3 = 1.53±0.47 of the CLF 
double power-law fits using Eq. 3, also plotted and labelled 
in Fig. 4 as Q2 and 03. 

We also performed a Schechter fit to the CLF of the 
target with the form of: 

cf>(M)dM = con x X as+1 e- x dM (5) 

where 

X = 10 0-4(M* -M) (6) 

which resulted in a characteristic magnitude of M£- « 
— 17.4 mag and a faint end slope as = 1.56 ± 0.19. Note 
that the values of a a and as are similar for the two types 
of fit. 

A bent LF for IRAS 18293-3413 fits better than a single 
power-law slope at the ~ 1.5a level. Though the difference 
is not very large, it is interesting, especially given that this 
galaxy has the largest statistics of SSCs in our present sam- 
ple. We will discuss the case further in Section 5.2.2. 



5 DISCUSSION 

5.1 The effect of blending on the LFs 

The distances of our sample galaxies, 45 < Di,(Mpc) < 100 
and one at 200 Mpc, means that any individual SSC cannot 
be resolved. Our spatial resolution of ~ 0. 1 " corresponds to a 
physical size of ~ 20 to 40 pc depending on the distance, and 
nearly 100 pc in the case of IRAS 19115-2124, and our pho- 
tometric apertures are of comparable size. Given that the 
effective radii of SSCs from the literature are in the range of 
3 to 5 pc (e.g. Whitmore et al. 1999), such aperture sizes may 
contain more than a single SSC candidate. Hence, blending 
of SSCs and complexes of SSCs will most probably contami- 
nate our SSC counts despite the use of AO imaging. Is it then 
reasonable to even refer to these as SSCs, or should we rather 
talk, for example, about "knots of star-formation" (Miralles- 
Caballero et al. 2011)? To address this we next attempt to 
estimate how much blending and crowding affect our anal- 
ysis. We first perform a simulation estimating the effect of 
blending on the LF slopes, then examine what happens to 
photometry of SSCs of a nearby system when it is moved to 
a larger distance. We also investigate the relation between 
the SSC surface density and confusion limits. In the end we 
conclude that blending effects are not significant within the 
SSC luminosity range considered, and in our case the term 
"SSC candidate" is a perfectly reasonable one for sources in 
our photometric catalogues. 

5.1.1 MC simulation of blending in LFs 

We performed a MC blending simulation to quantify the ef- 
fects of crowding on the values of the power-law index a. 
A random population of N artificial sources was created 
within the same magnitude range as our observational data 
(— 20 < Mk <■ —12) drawn from a LF with an initial index 
ainit = 2. We randomly selected two artificial sources from 
this population then blended them together. This process 
was repeated x times until the original population had x 
blended sources in total, corresponding to a blending rate 



y — x/N%. The blending was done step-by-step, so that 
a "new" blended source entered back into the catalogue as 
a new source and may be randomly blended with a third 
source, etc. A new power-law index, a new was determined 
from a fit to the blended source luminosity distribution. The 
difference in the slopes, Aa = ai nit — a new was then de- 
termined. We ran the simulation with y ranging from 5 to 
60 per cent in steps of 5 %. A particular blending realisation 
was repeated 1000 times for each blending rate. We also ran 
the MC simulation for ainit = {1-5, 1.8, 2.3}. We found that 
blends of 3 or more sources are rare below ~ 20 % blending 
rates. 

Figure 6, left panel, plots Aa against the blending rate 
with different values of ai„u. The sense is always that of 
flattening of the slope. A blending rate of 10 % for example 
would lead to a deviation of Aa ~ 0.14 from ainit = 2. The 
figure indicates that to get significant flattening of a by 0.3 
or more, one would need blending rates of ~ 25 % or higher 
for an initial slope of 2, or nearly 40 % for the slope to flat- 
ten from 2 to 1.5. The blending rates needed are higher for 
a steeper initial slope and lower for a flatter initial slope as 
expected. The blending rate is obviously related to the sur- 
face density of objects in the field, and we will return to the 
likely real blending rates of our targets in Section 5.1.4. 

Figure 6, middle panel, plots the change in magnitude 
of the brightest object; with more blending there is a greater 
chance that the brightest object is in fact a blend. Again we 
see that significant changes, Amag > 0.5, in the magnitude 
require blending rates of 40 % or higher. 

The simulation does not take into account varying de- 
grees of blending, nor changes in completeness limit. Most 
significantly, however, the simulation implicitly assumes a 
random surface distribution of targets, whereas SSCs are 
clearly clustered in (most) galaxies. While the simulation 
gives a good feel for the expected levels of changes in LF 
slopes due to blending, a more realistic estimate should be 
based on real data. 

5.1.2 SSCs in a redshifted Antennae 

The Antennae system, at _D_L~20Mpc, is a popular labo- 
ratory for SSC studies (e.g. Whitmore et al. 1999; Zhang 
et al. 2001), and we use its well-studied star cluster pop- 
ulation for a blending estimate. In particular, we retrieved 
//ST/WFC3 images (PI: Whitmore) in the UVIS/F814W 
(/-band) and IR/F160W (//-band) filters from the Hub- 
ble Legacy Archives. Adopting a distance of 22 Mpc the 
0.04" pixels in the /-band translate to a 4.3 pc physical size 
and the PSF to ~ 8 pc. In the //-band dataset the pixel size 
is equivalent to 9.6 pc and the PSF to ~ 16 pc. 

We first convolved and rebinned the /-band Anten- 
nae image to correspond to being 4 times further away. At 
88 Mpc it is well within the distance range of our sample (Ta- 
ble 1), and the PSF resolution element of 32 pc is also simi- 
lar to that in our NACO and ALTAIR/NIRI images. SSCs 
were detected both from the original and convolved images. 
We do not take into account the change in detection limits, 
so as to test only the effects of blending. The photometry 
and LF constructions were done similarly to our AO data 
described above, except that we do not attempt complete- 
ness corrections here as they are irrelevant to the main goal. 
Following other HST studies of the Antennae (e.g. Whit- 
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Figure 6. Results from our MC blending simulation. Left: the difference in the indices Act = cti n it — ct new plotted against the blending 
rate considering different values of cti„it- The figure indicates for instance that a blending rate of 8 % would lead to a deviation of ~ O.lf in 
the initial value of the index aiinti = 2. Middle: The same simulations, but plotting the magnitude difference of the brightest source. Note 
that intrinsic clustering characteristics are not taken into account in this simple simulation. Right: Surface density expressed as "beams 
per source" vs. blending rate; see Section 5.1.3 for details. A figure of 20 "beams per source" is traditional confusion limit rule-of-thumb, 
and it is seen that at this surface density the blending rate is ~ 10%, resulting in Ace ~ 0.1 to 0.2 depending on the initial LF slope, 
according to the left panel. 
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Figure 7. LFs of SSCs derived from HST images of the Antennae galaxies. Left: The black histogram shows the LF from the original 
/-band (F814W) image, and the blue after convolving and rebinning the image to correspond to a 4 times larger distance (w 84Mpc). 
The best-fit power-law slopes are indicated. Right: The same test for the H-band (F160W) image, but now the comparison is made so 
that the convolved image matches both the PSF size and the pixel size of our IRAS 18293-3413 NACO image. In both tests only mild 
flattening of the LF slope is observed. 



more et al. 2010) an aperture radius of 0.1" was used on the 
original (unbinned) image, corresponding to 10.8 pc. Nom- 
inal photometric zero-points and aperture corrections were 
taken from the WFC3 on-line manuals. The resulting LF is 
shown as the black points in Fig. 7, left panel, and a single 
power-law slope of a = 2.22 was found; the fit is performed 
down to Mi ~ —8.5 mag where incompleteness appears to 
be setting in. This value is in fact very consistent with other 
studies though no completeness or stellar contamination cor- 
rections were made in our test, for example Whitmore et al. 
(2010) found a = 2.26 in the V-band. 

The source extractions and photometry were done in 
the same way for the "redshifted" convolved image. An ad- 
ditional aperture correction is needed however since the PSF 
and pixel characteristics change, and the new correction was 
simply determined by matching the extracted magnitudes 
of a handful of foreground stars in the different images and 
apertures. The LF shown in blue in Fig. 7 results from ex- 



tractions from the convolved image when using a small 0.75 
pixel aperture, corresponding to 13 pc radius. As expected, 
the completeness limit is some 2 magnitudes brighter due 
to blending only, while the LF slope becomes only slightly 
flatter at a — 2.09. The aperture size used does not have an 
effect on the slope, but we will return to this aspect in more 
detail below. 

We also examined what happens to the 50 and 200 
brightest original SSCs in the convolved image. Of the 50 
brightest ones 7 were not detected as individual objects af- 
ter the "redshifting"(14% blending rate) while 53 of the 200 
brightest ones were not detected (27% blending). In the MC 
simulation of Sec. 5.1.1 these blending rates would have re- 
sulted in a flattening of the slope by Act « 0.20 and w 0.37, 
respectively. These numbers are somewhat larger than the 
measured flattening in the range 0.1 to 0.15 of the convolved 
image LF slope in the test above (Fig. 7). The recovered 
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magnitudes of the brightest 50 SSCs are a mere ~ 0.1 mag 
brighter than in the original image. 

In case clustering properties of SSCs would be differ- 
ent in the optical compared to NIR, we repeated a similar 
test using the //-band image of the Antennae. This time we 
modified the original HST image to match both the physical 
pixel scale and PSF size of the NACO IRAS 18293-3413 im- 
age at Dl ~ 75Mpc. Rebinning was not necessary because 
of the smaller pixels of the NACO instrument, while a fac- 
tor of 2.5 widening of the PSF was performed. Magnitudes 
were measured in 20 pc radii in the original image, matching 
the aperture used in the Whitmore et al. (e.g. 2010) NIC- 
MOS data, and 30 pc radii in the redshifted done 
for the NACO IRAS 18293-3413 images. The resulting LFs 
are shown in Fig. 7, right panel, indicating a flattening of 
the slope by 0.15, very similar to what was found for the 
/-band HST image. The magnitudes do not change signifi- 
cantly this time either, the average difference being less than 
0.1 mag for the 200 brightest SSCs. 

To understand and differentiate between the effects of 
spatial resolution and photometric aperture used with a 
given resolution, we redid the tests above with different con- 
volutions using PSFs ranging from 10 to 100 pc, as well as 
using numerous aperture sizes in the same range. Some re- 
sults become clear. First of all, the photometric aperture 
used at a given resolution does not change the slope signifi- 
cantly. Variations of ~ 0.1 were typically found for the value 
of q. Secondly, the LF does however shift to brighter mag- 
nitudes as the apertures grow. To recover as closely as pos- 
sible the intrinsic SSC counts the smallest possible aperture 
should be used, assuming a reliable aperture correction can 
still be determined. With the largest tested apertures the 
bright SSCs brighten by nearly a magnitude, while using 
apertures smaller than about 20 — 30 pc radius, the bright- 
ening stays within typical photometric errors of ~ 0.2 mag. 
Thirdly, when analysing the effects of different convolutions 
we found that the resulting LF slope remains relatively un- 
changed, staying within Aa = 0.2 of the original slopes 
of 2.15 and 2.22 (H and /-band respectively) until sizes 
of about 40 pc are reached, after which the slopes flatten 
rapidly reaching a ~ 1.5 at the poorest resolutions of a 
100 pc physical size. 

Figure 8 shows an example of a small region within the 
Antennae in both the original and convolved images of the 
HST /-band which illustrates these effects. While numer- 
ous faint SSCs disappear or are blended into bright SSCs, 
the latter are generally recovered in the convolved image 
with close to their proper magnitudes, unless located in very 
crowded regions. For example the more isolated bright SSCs 
on the edges of the image are recovered within ~ 0.1 mag 
of the original magnitudes, while the three bright SSCs in 
the middle of the central conglomeration have brightened 
by about 0.4 mag. The largest changes to the photometry 
are in cases where SSCs of similar brightness happen to lie 
near each other and blend into a single significantly brighter 
object, such as the bright SSC at bottom right which has 
blended into one source 1.1 mag brighter than either of the 
originals. How severe these effects are naturally depends 
upon on the overall surface density of objects, which is dis- 
cussed next. 
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Table 4. Surface densities, in units of kpc -2 , which would re- 
sult in confusion-limited observations of SSCs, when the limit is 
defined as 20 "beams per source". For small distances/PSFs the 
assumed 10 pc physical size of SSCs is resolved and the value cor- 
responds to this confusion limit; the value corresponding to the 
actual resolution element is given in parentheses. 

5.1.3 The effects of confusion limits 

Historically, blending and crowding as discussed above are 
referred to as "confusion" especially in the radio and far-IR 
studies, and in particular "bright source confusion". When 
a confusion limit is reached depends, in addition to the 
shape and size of the spatial resolution element, on the 
slope of the source counts. With very steep slopes confu- 
sion noise dominates, i.e. the many undetected sources at 
and below a detection limit give rise to an effective noise- 
level (see e.g. Vaisanen et al. 2001 and references therein 
for full discussions) that deeper observations cannot pene- 
trate. With much flatter slopes the bright sources lying too 
close together for satisfactory extraction at a given spatial 
resolution tend to dominate. The SSC source count slopes 
are closer to this latter regime, but it is prudent to search 
for this surface density-related confusion limit also in cases 
where the images are not truly confusion-limited yet. 

A 20 to 40 "beams per source" confusion limit is often 
used as a rule-of-thumb. For example, for a resolution ele- 
ment FWHM = 0.1", setting the beam size as 

FWHM 2 

then 20 beams per source corresponds to having one SSC per 
every 0.47" x 0.47" region. Hence, having over 500 sources 
within the area of IRAS 18293-3413 (see Fig. 1) would mean 
reaching the limit, while we detected ~ 200. However, in the 
core regions the surface density is definitely approaching the 
confusion limit. 

It is worthwhile to define a general confusion limit in 
terms of physical sizes for SSCs. Taking SSCs to be of size 
10 pc, the 20 beams-per-source criterion corresponds to a 
SSC density of 440 kpc -2 . If SSCs are seen more densely 
packed than this then better resolution is not likely to 
help extend SSC detection. Over the approximately 9 x 
9 kpc 2 size of the Antennae, for example, this would mean 
some 40000 sources; extrapolating the counts we extracted 
from the HST image, this level would be reached around 
Mi ~ —6.5 mag. As seen above, in individual regions of the 
galaxy system, the confusion limit must be reached at much 
brighter magnitudes since the SSCs are clustered. This is 
reflected in the fact that significant incompleteness starts 
appearing already at Mi ~ —8 mag in our test (Fig. 7, left 
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Figure 8. A 10" by 8"(1.1 by 0.9 kpc) region in the HST /-band image within the Antennae (region E in Whitmore et al. 2010). Left: 
The original image with ~ 160 SSC candidates detected. Right: The image after taking the galaxy four times further away. Only ~ 20 
SSCs are now detected due to blending, though the brighter population remains relatively unchanged. 



panel), //the HST image was truly confusion-limited, and 
assuming a LF with a — 2, the completeness limit would 
be expected to brighten by 3 magnitudes after the factor 4 
convolution, whereas a change of 2 mag was observed. 

We list in Table 4 physical surface densities per kpc 2 
at which the confusion limit is likely to be reached at a 
given distance and spatial resolution, given the definitions 
above. With the typical PSF of ~ 0.1", the distances of 
our sample, and the number of detected SSCs, only the 
core regions of IRAS 18293-3413 come anywhere close to 
these confusion limits, as well as IRAS 19115-2124 due to 
its distance. This is not to say that individual SSCs will not 
blend of course, especially in clustered star-forming regions 
of the galaxies as shown in the tests based on the Anten- 
nae. For example, the surface density of detected objects in 
the left panel of Figure 8 is ~ 160kpc~ 2 . At a distance of 
~ 20Mpc and for a resolution element of 0.075 "this is not 
quite yet confusion-limited according to Table 4, though it is 
approaching it. This is also seen from the corresponding 78 
beams-per-source surface density. However, within the clus- 
tered sub-region the equivalent values would be ~ 500 kpc -2 
and ~ 25 beams-per-source, respectively, i.e. that region is 
confusion-limited. 

To see how the surface densities quantitatively relate to 
the LF slope changes, we ran another MC simulation adding 
increasing amounts of randomly distributed equally bright 
sources in an otherwise empty frame and extracting them 
with SExtractor using typical parameters. This was done 
hundreds of times at several surface densities, and the frac- 
tion of unrecovered sources is simply the blending rate at 
each surface density. The change in LF slope with a given 
intrinsic LF shape was already simulated in Section 5.1.1 as 
a function of this blending rate. The right panel of Fig. 6 
connects the two by plotting the "beams per source" surface 
density against the blending rate. The black curve is the in- 
trinsic surface density of the simulation, and the red curve 
is calculated from the extracted surface density. An ob- 
served surface density of objects can be used to get an ex- 
pected blending rate using the red curve (for completeness- 
corrected counts the black curve is more appropriate) and 



this can then be converted to a likely Aq value of the LF. 
While the vertical displacement of the surface density vs. 
blending rate curves will depend on the source detection al- 
gorithm and clustering of objects, our tests show that the 
relations do serve as a realistic approximation of the quan- 
titative effects involved. 



5.1.4 Effect of blending in our sample 

The average surface density of detected SSC candidates in 
our target galaxies ranges from a low of 0.3 kpc -2 in IC 694, 
or 168 beams per source, to a high of 0.8 kpc -2 in the case 
of IRAS 18293-3413, or 87 beams-per-source. According to 
Fig. 6, the corresponding blending rates are well below 5 %, 
meaning that Aa < 0.1 for any LF shape, most likely < 0.05 
for most of the targets. Because of the uncertainties involved, 
such as the difficulties of properly accounting for clustering, 
we prefer not to attach Aq values for each individual target, 
but rather stress that these must be less than 0.1. However, 
in our most extreme case of the central ~ 2 x 2 kpc region of 
IRAS 18293-3413, the SSC surface density equates to 20-30 
beams per source, and Aa ~ 0.15 would be expected in that 
region. 

As a further check, Fig. 9 shows separately the LFs for 
SSCs out to distances Dl ^ 60Mpc, and for Dl > 60Mpc 
(excluding IRAS 19115-2124 since it is so much further 
away). If blending was a serious issue in our sample, the 
more distant subsample would be expected to show a flatter 
slope. The single power-law indices of the LFs below and 
above 60Mpc are a = 1.98 ± 0.11 and a = 1.89 ± 0.23 re- 
spectively, i.e. there is no significant difference in the LF 
slope, nor in the numbers of detected SSCs, with distance. 

In summary, we are confident that in the luminosity 
range of interest, Mk < —14 mag, any flattening of LF 
slopes due to blending is below Aq ~ 0.1. With the possible 
exception of IRAS 19115-2124 where the spatial resolution 
is just a physical size of 90 pc, the luminosities of detected 
if-band point sources will be dominated by a single bright 
SSC rather than whole knots of star formation. In the case of 
Antennae-like SSC populations, blending and crowding do 
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Figure 9. SSC LFs of two subsamples generated from our observational data but segregated by distance. The dataset from IRAS 19115- 
2124 is not included to avoid bias in the analysis. Left: SSC LF of the closer subsample where the targets have distances ^60Mpc. 
Right: SSC LF of the more distant targets, excluding IRAS 19115-2124. The values of the slopes appear to be consistent within the 
uncertainties. 



flatten the LF slope, but significantly so only at resolutions 
poorer than a 40 pc physical size. The photometric apertures 
used should be as small as possible to recover intrinsic lu- 
minosities, though the aperture does not have a significant 
effect on the LF shape. Assuming a 10 pc scale for SSCs, 
the confusion limit is reached at a surface density of 440 
SSCs per kpc 2 , or less when clustering and larger resolution 
elements are a factor (Table 4). 

5.2 Interpreting the values of a 

The results from our fitting procedure in § 4 show that a 
single power-law distribution is a reasonable approximation 
for the young stellar cluster LFs in our sample. A single 
power-law SSC LF has been reported by many theoreti- 
cal and observational studies (e.g. Elmegreen & Efremov 
1997; Lada & Lada 2003; Gieles et al. 2006a), and such 
behaviour may be linked to the shape and characteristics 
of the physically more interesting underlying (initial) mass 
function and/or the evolution and disruption of the SSC 
population. As listed in Table 3, the fitted single power- 
law indices range from a ~ 1.5 to 2.4, with average sin- 
gle slope values of a con = 1.87 ± 0.30 for a constant bin 
size, and a var — 1.93 ± 0.23 for a variable bin size. Val- 
ues of a reported in the literature vary widely from 1.7 to 
2.4 (e.g. Elmegreen & Efremov 1997, Whitmore et al. 1999, 
Elmegreen et al. 2002), with slopes in normal spirals often 
in the upper part of this range. Our a values tend to be 
in the less steep part of the range, and in at least one case 
(IRAS 18293-3413) either a broken power-law or a Schechter 
function would yield a better fit. These two issues are dis- 
cussed next. 



5.2.1 LF slopes: a difference between LIRGs and normal 



Improper incompleteness corrections can produce artificially 
flattened slopes. We note however that the slopes are fitted 
to a range where at most 20 % of SSCs are missed, and where 



the corrections are still reliable. The most significant change 
is likely to be caused by blending of SSCs but as discussed 
in the previous section, this is expected to be Aa ~ 0.1 
at most for our targets, and likely less. Also, we did not 
see any significant change in the a values with distance to 
the galaxies, which would have been expected if significant 
blending was present (Fig. 9). Therefore, we conclude that 
the intrinsic SSC LF slopes, the combined and averaged LFs 
having a ~ 1.82 ±0.25 and 1.87 ±0.30, respectively, are not 
affected by these observational effects, with any systematic 
effects being outweighed by the statistical uncertainties. 

Normal spiral galaxies tend to have power-law indices in 
the range of a ~ 2 to 2.4. Calculating the average from SSC 
LF studies which include several spirals each (Larsen 2002; 
Gieles et al. 2006a; Mora et al. 2009), we get a « 2.17±0.24 
from 10 spiral galaxies. If we include other individual galaxy 
studies (Schweizer et al. 1996; Buta et al. 1999; Elmegreen 
et al. 1999, 2002; Cantiello et al. 2009; Pellerin et al. 2010; 
Whitmore et al. 2010), where some targets are starbursts, 
though not yet LIRGs, we get an average of a « 2.18 ± 0.19 
from a total of 19 galaxies. Note that the values come from 
different filters, mostly from V- and /-bands; if several bands 
were used, we chose the reddest band. As discussed earlier 
the well-studied interacting Antennae system also has a ~ 
2.2 (Whitmore et al. 2010). 

While most of the previous SSC LF determinations in 
the literature, including the Antennae, have been in actively 
star-forming galaxies, the vast majority of them have milder 
SFRs than LIRGs. Our average power-law slope for LIRGs, 
a ~ 1.9, appears flatter than the slope in normal galaxies. 
The large sample of LIRGs studied in Vavilkin (2011) shows 
an average a ~ 1.8, though blending effects there might 
still need to be explored. Similarly, Miralles-Caballero et al. 
(2011) find a ~ 1.9 for their (U)LIRGs at < 100 Mpc, and 
even natter values at certain interaction stages of LIRGs. 
Moreover, Adamo et al. (2010, 2011b) recently found flatter 
power-law slopes while probing the star cluster properties in 
Haro 11 and Mrk 930. Both of the targets are blue compact 
dwarfs with intense star formation, and the former can also 
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Figure 10. SSC LFs of two subsamples generated from our observational data, excluding IRAS 19115-2124, but broken by the average 
SFR. Left: SSC LF of the subsample where the targets have SFR sg61MQyr _1 . Right: SSC LF of the other set of targets with high 
SFRs. A power-law fit until the 80 % completeness level generates a quite similar value of a in both cases. 



be classified as a LIRG. Is it the case therefore that LFs are 
systematically flatter in extreme SF cases? 

To see if there is a trend with SFR within our sample, 
Fig. 10 shows LFs separately for SSCs with host galaxy SFRs 
less than, or greater than the average SFR~ 61 MQyr~ l , 
again with IRAS 19115-2124 excluded. The values of the 
slopes are a = 1.89 ± 0.16 for the subsample below the 
average and a = 1.88 ±0.19 for those above. Thus, while we 
do not find a correlation between the LF slope and the SFR 
within our present small sample, our results together with 
other recent studies suggest that there is a real difference 
between the SSC LF slopes of LIRGs and those of more qui- 
escent galaxies. This needs to be verified with larger samples 
for better statistics, and with careful blending analysis since 
larger samples will by necessity involve more distant targets. 

Apart from effects related to the observations them- 
selves, perhaps the most fundamental cause for flattened 
LFs, or breaks/bends in the LF for that matter, would come 
from mass- and/or age-dependent cluster disruption, as well 
as from differences in the cluster formation with different 
environments. Young star clusters are most vulnerable to 
disruption, leading to variation of the integrated LF as time 
passes. Since LFs are the integrated sum of the distribu- 
tions of individual initial LFs of SSCs of different ages and 
masses, any selective (e.g. mass-dependent) disruption of 
them would be seen as changes in the (integral) total SSC 
LFs (e.g. Mengel et al. 2005; Gieles et al. 2006a; de Grijs 
& Parmentier 2007). Kruijssen et al. (2012) argue that clus- 
ter formation efficiencies decrease in galaxies with higher 
SFRs due to increased tidal disruption in the extreme envi- 
ronments ("cruel cradle effect"), which might lead to flatter 
LF slopes if lower mass SSC formation is more affected. Al- 
ternatively, high SFRs may favour the formation of more 
massive GMCs (e.g. Wei et al. 2012). However it is impossi- 
ble to conclude whether disruption is the main cause of the 
flatter slopes until we have analysed data from other filters 
to be able to derive LFs and MFs in different mass and age 
regimes. 



5.2.2 IRAS 18293-3413: What could a broken power-law 
reflect? 

In the case of IRAS 18293-3413 fitting two independent 
slopes produces better Chi Square values (Xred = 0.57) 
than a single power-law. The break appears to be at Mk ~ 
— 16.5 mag and persists regardless of how the data is binned. 
We note that the distribution can also be well fit with 
a Schechter function, giving a characteristic magnitude of 
M K —17.4 mag (Xred = 0.64). In general, any bends or 
cutoffs found close to the detection or completeness limits 
should be viewed with caution, but in this case the break is 
more than a magnitude above the 80% completeness limit. 

Bends have been argued to result from physical trunca- 
tions of the SSC high-mass end of the SSC mass functions 
or from disruptions of the less-massive SSC population (see 
e.g. Gieles et al. 2006a, Haas et al. 2008). The former results 
in a steepening of the bright end of the LFs if the truncation 
is a physical one and not merely a result of a size-of-sample 
effect; the latter affects the faint end of the LF by making it 
flatter. Thus both can result in bends in the shape of the LF, 
making for better fits with broken power-laws or Schechter 
functions. The absolute magnitude of the break point would 
increase with the maximum mass of the SSC population, 
though the exact relation is a complex function of the dis- 
ruption models. The break is expected to occur at brighter 
magnitudes at redder wavelengths because of the expected 
red colour of the aged SSCs. Whether or not these bends 
and truncations are in fact seen in LFs or mass functions is 
still a matter of intense debate (see e.g. Gieles et al. 2006a, 
Larsen 2009, Whitmore et al. 2010, Fall & Chandar 2012 
and references therein). 

Bends in SSC LFs have typically been found in the 
range of M ~ —9 to —11 mag in VRI filters (e.g. Whitmore 
et al. 1999; Gieles et al. 2006a,b; Santiago-Cortes et al. 2010; 
Adamo et al. 2010, 2011a). Assuming typical un-extincted 
(and age dependent) optical-to-if-band colour indices of 2 
to 3 mag, one would expect a corresponding bend to appear 

at M K ' 12 to -13 mag. Mengel et al. (2005) find a bend 

at a significantly brighter level Mk ~ —15.5 mag in the An- 
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Figure 11. Starburst99 model of 2 X 10 5 Mq SSP cluster with 
a Kroupa IMF. The three curves show the absolute magnitude 
of the star cluster in B-(blue), /-(green) and K-(red) bands. The 
inset highlights that the X-band luminosity of star clusters is 
likely to peak close to age lOMyr while they peak somewhat 
earlier in the optical. 

tennae (but see Fall et al. 2009) not unlike our result for 
IRAS 18293-3413, suggesting significant extinction effects. 
The bend-point of a double power-law at Mk ~ —16.5 mag, 
or the Schechter characteristic magnitude at ~ —17 .4 
mag, would correspond to a minimum characteristic mass of 
M* r; 1-2 x 10 6 Mq , or perhaps double this mass range if 
we adopt an average mass-to-light ratio over the first 30 Myr 
of age (see next Section). These masses are higher by fac- 
tors of 5 to 10 than the characteristic mass found for SSCs 
in normal spirals (Gieles et al. 2006a; Larsen 2009) while 
they are very similar to those suggested for galaxies with 
much higher SFRs, such as LIRGs (see e.g. Bastian 2008). 
However until there is clearer evidence for real bends in SSC 
LFs from a larger dataset it is premature to read too much 
into these masses. 



5.3 Mass range approximations of the I-f-band 
SSCs 

An interesting characteristic of SSCs observed at NIR wave- 
lenghts is a relatively narrow time-frame at ages of ~ 10 Myr 
when the high mass stars in the SSCs enter the RSG phase 
making them very NIR-luminous and "suddenly" red in the 
optical-to-NIR colours compared to earlier blue stages. This 
is seen in Fig. 11 which plots as an example the evolution 
of cluster brightness with time in the BIK filters, derived 
from a Starburst99 (Leitherer et al. 1999) model assuming 
an instantaneous SF with a fixed mass of 2 x 10 5 Mq, and a 
Kroupa IMF. Hence, a lower limit on SSC masses detected 
in this work can be obtained by assuming the mass-to-light 
ratio at that ~ 10 Myr age. In many of our galaxies the most 
massive detected SSC has Mk ~ —18 mag corresponding 
to a mass of 4 x 10 6 Mq at 10 Myr of age or more than 
2 x 10 7 Mq at 30 Myr; clearly we are sampling very massive 
clusters in LIRGs. 



SSCs with masses in excess of 10 7 Mq have been found 
before (e.g. Bastian et al. 2006; Portegies Zwart et al. 2010). 
The faintest-detected SSCs in our sample have Mk ~ — 13 
corresponding to masses of « 2 x 1O 4 M , while the pho- 
tometric 80 % completeness limits in different targets corre- 
spond to lower limits of 1 — 2 x 10 5 Mq . We note that with 
such massive clusters, the luminosities of individual bright 
RSG stars are negligible compared to the total integrated 
flux of the star cluster. Therefore, stochastic effects should 
not introduce significant scatter in the mass-to-light ratios 
and inferred ages of our clusters (Fouesneau et al. 2012). 



6 SUMMARY 

In this work we studied the characteristics of massive star 
clusters in the extreme environments of local interacting 
LIRGs. Observations were performed in the if-band filter 
using two different NIR AO instruments with pixel scales 
of ~0.02"pix _1 and a FWHM~0.1" for point sources. 
The galaxy sample consists of LIRGs in the redshift range 
0.01 < z< 0.05. This study is opening new territory in SSC 
studies since both the distance range and host galaxy SFR 
range are higher compared to most previous SSC studies. 
Specifically, we derived the K-h&nd luminosity functions of 
SSC candidates in our targets, and because of the distances 
involved, we also carefully evaluated the effect of blending 
on the power-law index a of the LFs. Extensive photomet- 
ric completeness simulations were done, as well as checking 
the effects of sample binning and foreground contamination, 
which turned out to not be significant. The main results from 
this work can be summarised as: 

(i) The SSC luminosity function is probed at high com- 
pletion down to Mk ~ —14 or —15 mag in our sample. In 
this range all the LFs are reasonably well fitted by a single 
power-law, though in the case of IRAS 18293-3413 a double 
power-law or a Schechter function is a better approxima- 
tion. The values of the best-fit slopes vary with a wide range 
from a = 1.5 to 2.4 with the average value at a w 1.9, and 
the combined SSC LF, excluding the most distant target, 
at a w 1.8. The slopes appear slightly flatter than those in 
normal spirals which typically have a ~ 2.2. The sample 
is still too small to find definitive trends of the slope with 
SFR or interaction stage. One or more of age, extinction, 
and mass-dependent cluster disruption effects can all lead 
to small a values at the faint end, but cannot be unambigu- 
ously separated from the present dataset alone. 

(ii) We carefully examined the possibility of blending of 
SSCs in our target LIRGs. Though blending does happen 
with the typical resolutions of 30 to 40 pc physical sizes of 
our sample, we showed that it is not enough to change the 
LF slopes by more than Aa w 0.05 — 0.1, nor change drasti- 
cally the measured luminosities of SSCs. Hence we conclude 
that out to ~ 100 Mpc it is quite possible to accurately mea- 
sure SSC properties with 8-m class telescope adaptive op- 
tics and HST imaging. In addition to deriving some general 
blending/confusion properties, we found that the photomet- 
ric apertures used do not affect the LF slope, but small aper- 
tures are necessary to recover the luminosities as correctly 
as possible. Worsening spatial resolution tends to flatten the 
measured LF slopes through increasing blending. However, 
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the effect becomes pronounced only when close to the con- 
fusion limit. In the case of SSCs distributed as in the An- 
tennae, we determined that LF slopes at Mr < —12 are not 
reliable if the spatial resolution corresponds to a physical 
size larger than ~ 50 pc. 

All our findings are based on observations with a single 
filter and thus the estimation of ages of the star cluster can- 
didates was beyond the scope of this work. However, we are 
in the process of estimating the physical characteristics of 
the selected star cluster candidates with the help of archival 
HST data and by using J and H imaging from Gemini to be 
able to constrain the SSC evolution models more robustly. 
In addition, with a larger sample of southern LIRGs and 
starbursts with the VLT/NACO currently available, and the 
next generation of multi-conjugate AO systems promising to 
deliver much more stable PSFs across a larger field of view 
we expect to grow both the number and quality of SSC LFs 
to probe correlations with LIRG host galaxy characteristics 
and environments. 
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